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Abstract 

A functional method to achieve the summation of all Feynman graphs relevant to a par- 
ticular Field Theory process is suggested, and applied to QED, demonstrating manifestly 
gauge invariant calculations of the dressed photon propagator in approximations of increas- 
ing complexity. These lead in a natural way to the extraction of the leading logarithmic 
divergences of every perturbative order, and to a demonstration of the possible cancella- 
tion of all such divergences in the calculation of the (inverse of the) photon's wavefunction 
renormalization constant Z3. This analysis provides a qualitative understanding of why the 
measured value of the renormalized fine structure constant is, approximately, 1/137. 
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1 Introduction 



Generations of Physicists have grown to maturity with the aid of Feynman graphs for 
various perturbative problems; and when the subjects of their theoretical researches have 
required a non-perturbative approach, they have attempted to sum sub-sets of Feynman 
graphs in order to build some semblance of a non-perturbative approach. 

The only trouble with this approach is that it is physically and mentally impossible 
to sum all Feynman graphs needed for a particular amplitude. Consider, for example, a 
simple process of electron-electron scattering in QED for which one contemplates summing 
the contributions of all relevant Feynman graphs; and for simplicity let us neglect those 
" vacuum " graphs with closed electron loops, as well as self-energy graphs for each elec- 
tron, and consider only graphs built from the exchange of virtual photons between the 
two scattering electrons. As each order of perturbation theory is increased, a new class of 
graphs is generated, none of whose members can be obtained by iteration of the previous 
graphs; and as the perturbative order increases, there will be, eventually, an infinite num- 
ber of graphs in each class. Bethe-Salpeter equations, whose kernel contains elements of a 
finite number of classes of graphs, are a so-called " non-perturbative " attempt at bulling 
a solution which contains all powers of the coupling; but it cannot yield the correct answer 
because it has missed all the classes of graphs not contained in the kernel. And there is no 
known non-perturbative method of calculating that kernel. 

Is there an analytic method of summing all relevant Feynman graphs ? The answer 
is positive, but it is a functional method, made possible by two distinct developments, 
both of which have been known for many years : the Schwinger-Symanzik [Ij functional 
expression for the Generating Functional (GF) of the Quantum Field Theory (QFT) in 
question; and the functional representations first given by Fradkin [2] of the two functionals 
contained in that GF. The functional operations which yield all possible n-point functions 
of the theory are equivalent to a Gaussian-weighted functional integration over repeated 
powers of those two Fradkin functionals, depending on the process contemplated and on the 
degree of precision in which one is interested; an alternate and somewhat simpler " linkage 
operator " formalism can be employed which performs the same functions. In either case, 
the functional operation is defined over appropriate fluctuations of a space-time field A(x), 
which are defined precisely by the Fradkin functional representations. 

All of perturbation theory, which in any order is identically equivalent to the sum of all 
Feynman graphs of that order, may be obtained by a simple, perturbative expansion of the 
GF in powers of the relevant coupling constant; in any finite order of expansion of the GF, 
the functional operations may be performed exactly. But can the needed, non-perturbative 
functional operations be performed ? The answer is again positive, because the Fradkin 
representations are themselves Gaussian in their A dependence, and an exact evaluation of 
a Gaussian weighted integration over Gaussian functional dependence can be performed. 
This effective summation over all Feynman graphs then produces an answer given in terms 
of the Fradkin functional representations of the two functionals alluded above; and one may 
well ask : Is this progress ? 

Again, the answer is positive, and the reason is that the Fradkin representations are 
Potential Theory constructs, with decent approximations in different physical situations. 
In QED, for example, which QFT we shall here use to illustrate these remarks, the two 
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functionals whose Fradkin representations are needed, are the electron Feynman (or causal) 
propagator Gc{x,y\A) in a specified vector potential A^{x) : 

[mo + 7-(5x - igoA{x))]Gc{x, y\A) = d^^^x - y) 

or, defining the operator Gc[A\ by : <x | Gc[A] \ y> = Gc{x^y\A) : 

Gc[A] = [mo + j-{d-igoA)r' 

and the " fermion determinant ", or " vacuum functional " e , where : 

L[A]=TTln[l-igoj-ASc], 5, = G,[0] 

This fermion determinant represents the sum of all, single fermion loop graphs which contain 
all possible ( even ) numbers of attached photon lines. 

The relativistic notation used throughout this article is the so-called "east coast metric " , 
with the scalar product defined by a-b = a-b—aobo, and the Dirac matrices such that = J^, 
7^ = 1 and { 7^, 7,, } = 25^^. 

At high energies, a very good approximation to Gc[^] is a Bloch-Nordsieck, or eikonal 
approximation [3], and its various corrections, which are also ( and even simpler ) Gaussians. 
Corresponding approximations may be defined for L[A]; but it is the special, exponenti- 
ated way in which L[A] appears which defines the nature of useful approximations; e.g., a 
functional cluster approximation. 

These statements will here be illustrated by successive calculations of higher order cor- 
rections to the photon propagator of QED, ending with a simplified, intuitive extraction of 
the sum of the most relevant, divergent terms in the expression for the inverse of the pho- 
ton's wave function renormalization constant. Gauge invariance in every order is manifest. 
In contrast to previous approximations which were unable to sum exact forms [3], we find 
a mechanism which has the possibility of producing a cancellation of the divergent pertur- 
bative logarithms, suggesting that Z^^ may actually be finite, generating a finite charge 
renormalization; in the process, we find arguments as to why the renormalized fine struc- 
ture constant is on the order of 1/137, independently of the mass of the fermion ( electron 
or proton ) possessing that electric charge. The latter calculation is intuitive, rather than 
rigorous; but it is a compelling indication of the power of the functional approach described 
here. 

To be more explicit, we briefiy review in Section 2 the Schwinger-Symanzik functional 
expression for the Generating Functional (GF) of QED and the convenient version of the 
Fradkin representation for the functional fermion determinant. Then, in the context of this 

functional approach, we give a simple derivation of the lowest order vacuum polarization 

~ (2) 

Kliu{k). That calculation is then generalized in Section 3 to display the sum of all such 
KfnJ^\k), corresponding to the exchange of all possible virtual photons across and on the 
same fermion lines of that simplest closed fermion loop graph; and in the process we discuss 
two sets of special cancellations which simplify all subsequent calculations immensely. ( In 
contrast to the usual Feynman graph calculations, the cancellation of certain divergences is 
apparent immediately, before the relevant functional integral ( FI ) is evaluated ). We set 
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up the old Jost-Luttinger [5] calculation of itju) (k) in this new and simpler formalism, and 
leave it as an exercise for the interested reader. 

Each of the Kju!^\k) so calculated displays one or more log divergences, as does the 
photon's Z^^. We introduce in Section 4 the Dominant Part ( DP ) Model, which extracts 
all possible log divergences from every order of the graphs defining 



:i.i) 



n=l 



and indicate how and why the result is still divergent. At this stage, perturbative renor- 
malization can be carried out in the conventional way. 

But the complete K^y{k) contains an infinite number of closed fermion loops, with all 
their radiative corrections, attached in all possible ways to the graphs which gave (1.1); and 
in Section 5 we introduce an Extended DP Model which extracts and sums the divergent 
parts of all graphs which contribute significantly to Z^^ . With the use of simplifying 
but reasonable approximations to the two resulting integrals, we then find that the Z^ 
so obtained is finite. The requirement that the ImZ^ = then defines an " eigenvalue 
condition " for oq = go/4:-jr; and together with the finiteness of Z^ , the possible value(s) 
of a()Z^{ao) define the renormalized a. Within the context of our seemingly reasonable 
approximations, it follows that ^3^^(00), and hence (oq), are independant of the fermion 
mass; and that the expected value of a should be close to 1/137. 

Finally, in the Summary of Section 6, we compare this analysis with previous, unsuc- 
cessful ones of four decades ago. We emphasize that the arguments and proofs given in 
this paper are intuitive and heuristic; they are " physicists proofs ". We believe that they 
provide a framework in which one can resolve the long standing puzzle of whether QED, 
and other gauge theories, can be considered as finite and true theories of Nature. 



2 The Generating Functional and Fradkin's representation 

We begin by stating the Schwinger-Symanzik Generating Functional for QED, defined 
in a covariant gauge : 



rQED 



1 



<0|5|0> 



exp 



igo d X- 



5ri{x)\ 6Jt^{x) J 5rj{x 
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(2.1) 



exp 



i d'^x d^y 7){x) Sc{x - y) r^{y) + - d'^x d'^y J^'ix) D,^^,{x - y) r{y) 



An alternate, equivalent, and somewhat more convenient representation is given by : 



1 



exp 



- ld^xd^yJ^'{x)D,^^,{x-y)r{y) 



<0|5|0> 

xe^-expf^/d^.Ar/(x)G.(x,,|A),(,) exp[L[^]] 
(x-y) r{y). 



(2.2) 



with A^{x) = Jd'^yDc^fj.u 
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The normalization Z^^^[0,0,0] = 1 is defined by : <0|S'|0> = e^^ exp[L[yl]] 
and e'^^ represents the " hnkage operator " : 



A=o' 



{x - y) 



The free photon propagator is given in the Feynman gauge by : Dcuuik) = —rr— — 

/c^ — ie 



m 



momentum space, and Z?c 



in configuration space. 



In a previous paper deahng with both photon and electron propagators [6], use was made 
of the Fradkin functional representation for Gc{x,y\A), in the calculation of the dressed 
electron propagator, but the calculation was performed in a " quenched approximation ", 
negelecting effects of the vacuum functional L[A\ : 



S',{x-y) = e'^AG,[x,y\A) 



<S> 



A=0 



e^^Gcix,y\A) 



A=0 



In this paper, we treat radiative corrections to the fully dressed photon propagator : 
Dc,^ll'ix -y) = Dc,f,uix -y)+ 11 Dc^^xix - u)Kxa{u - w)Dc^„^{w - y)du dw 



iK^u{x - y) 



6^L 



6Af,{x) SA^{y) 
5L 5L 



> 



+ 



_SA^{x)6A^{y) 6Af,{x) 6A^{y) 
in contrast to the quenched approximation of ref.[6] : 

6^L 



A=0 

/<S> 



A=0 



iK^y{x - y) 



5A^{x)5A^{y) 



A=0 



(2.3) 



(2.4) 



and find that our conclusions are strongly dependent upon the effects of exp[L[A]]. 

We here sketch the evaluation of the lowest order Kl^J (k) directly from the gauge invari- 
ant L[A] formalism, obtained from (2.4) without the linkage operations of that equation. 
The fermion determinant has an exact Fradkin representation of the form : 



L[A] 



1 



2 Jo s 
-igo I ds' Vf,{s') Aij,{x' 



2 ' 
-isiriQ g 



ds' 



<5(^)(/ ds'v{s')) jd'x' 



X e 



v) 



tr e 



rs rs' 

go / ds afj,uF^y{x - I v) 
; Jo Jo 



v=0 



Lo (2.5) 



where F^y = d^j,Ai, — di^A^, a^y denotes the gamma matrix combination \[^^-,'~^v\, Vf^{s') 
denotes the four velocity of the fermion with an instantaneous proper time s' and the 
constant Lq is such that L[0] = 0. A more convenient version of this representation is 
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obtained from the variable change u^{s') = Jq ds"v^{s"), which yields the normalized 
functional integral : 



L[A] = -\f'ie ---^ Jd^x' N jd[u] e \ J "^'"^^^^ S^^^ 

rs rs 

-igo / ds' u' Js') Af,{x' - u{s')) /go ds' a^,^ F^,„{x' - u{s'))\ 
xe J'i tr e ^0 ] - Lq (2.6) 



where N = d[u\ e2j , <si\h \s2> = -7: — S{si — 82)-^ — ^^'i h( 81,82) = 

J OSi O82 

inf(si, S2) = ^(si + S2- \si - S2I) = si 9{s2 - si) + 82 6(81 - 82). 

It will be most relevant to remind the reader of certain rigorous properties of L[A], 
which quantity in fact depends only upon F^^, and can be easily be written as such . 
This property immediately carries the consequence that currents induced in the vacuum, 
< Jfii^) > = ig ^L/SA^{x), are to satisfy charge conservation : jfj,{x) > = 0. In terms 
of the K^i, of (2.3) or (2.4), this means that d^K^^^ = di,K^y = 0, so that the K^i,{k) 
are expected to have the gauge invariant form {k^ky — k"^ 5 ^j,y)Ii{k'^) . The simplest, order 
a = 5^/47r, Feynman graph corresponding to a single closed fermion loop does not display 
this property; and in the past, special, ad hoc maneuvers were invented to restore gauge 
invariance. In the Pradkin representation for L[A\ used in this paper, gauge invariance to 
all orders is automatically satisfied. 

It is important to note that there are two restrictions on the u{s') variables, the first 
an implicit condition u(0) = 0, which arises from the definition of u{,s'); and the second 
condition, u(s) = 0, explicitly stated by the delta function of (2.6). For this lowest-order 
calculation, we replace mo by m. 

Performing the pair of functional derivatives on L[A] yields : 

xAgl j'dai [ds2 -<(siK(s2) - (S„8„ - V*^) i^lx - » + - ufe)) (2.7) 

after making use of the trace properties tr(l) = 4, and tr((7Q^(7^i/) = Sav^p^j. — ^ap^^u- 
Writing a Fourier representation of 5^^\x — y + u(si) — u{s2)), the second line of (2.7) 
becomes : 



fd'^k / dsi /ds2e^^-(^-^/) + ^^-(^(*i)-"(^2)) 
(27r)4y Jo Vo 

-u'^{8i)ul{82) + {k/^ky - S,^„k'^) 



(2.8) 



and one can see that gauge invariance is maintained by imagining the result of calculating 
djj,, or equivalently, of multiplying (2.8) by A;^ and summing over /x. The term coming from 
the sigma matrices obviously vanishes; and the u'{8i) dependence yields : 
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Jo JO OSi 

because tt(0) = u{s) = 0. At this early stage, it is then clear that the result of this 
computation must be gauge invariant. 

In order to avoid any confusion on taking the si^2 derivatives of the u' factors with 
the u{si) — u{s2) exponential terms, it will be useful to introduce the sources = 
kfj,[S{s' — si) — 6{s' — S2)], and also /^(s') = 6(s' — s), and there then follows : 



d d 



g —ism 



s 
5 



dsa dsb dgf,{sa) dgu{sb) 



N j d[u] e 2 
i fu.{f + g) 



u(2/i)"^u 



(2.9) 



The normalized FI over the u dependence is then immediate, and yields the exponential 
factor : 



so that the combinations 



e 

d d 



-i JJ\f + 9)Mf + 9) 



(2.10) 



generate : 



dsa dsb Sg^isa) Sgt,{sb) 

-4bM + knidisi - Sa) - e{s2 - Sa))\\pu + k^{e{si - Sb) - 9{s2 - Sfe))] 

-2iS,j,uS{sa - Sb) 
and in the limit as Sa ^ si, Sf, — )■ S2, produce : 

{-2i)5ij,,yS{si - S2) - ^IpixPu + ^k^ku + {p^lku + Pukfj,){9{si - S2) - ^)] 
With the exponential factor of (2.10) evaluated as : 

- isp^ - 2isi2p.k - iA;^|si2| , S12 = si - S2 



(2.11) 



(2.12) 



(2.13) 



the / d^p may be performed, leading to the cancellation of the term coming from the 
tr[o"a//0"/3i/] dependence, and the result : 



87r2 



3/7e ■ 2 1 



s 



2% 8^,- - 2i 5^, 5{si -S2) + A (1 
s s 



Si - S2\ 



I kuki/ Sui/k 



X exp 



Si - S2 



(2.14) 



Finally, with the aid of the relations : 

rs rs rl 

/ dsi / <is2/(|si - S2I) = 2s^ / dy{l-y)f{sy) with y 
Jo Jo Jo 



Si - S2 
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f'dy (1 - 2y) e -i^k^^ " = ['^V |- e "^^^'^(1 - v) = 
Jo sk^ Jo dy 

f'dy y{l - 2y) e - y) = _ f\ e -^^^^{1 - y) 

Jo sk^ L Jo 



one obtains the result : 

K^Kk) = {k,k,-5,uk^)^^'\k') , n(2)(fc2)=/ ^e-^^"^ n(2)(fc2,.) (2.15) 

Jo s 

in which the gauge symmetry has been preserved, and the subsequent / ds leads to a log 
divergence at its lower limit. With : 

n(2) {k\ s) = ^j\yil-y)e -^'^Ml " v) 

one has the simplest vacuum polarization result of a half-century ago [7J. Then, with the 
inverse of the photon's wave function renormalization constant given by : 

z^^ = 1 + n(2)(o) 

the renormalized, to order a, vacuum polarization is given by : 

n(2)(^2)_n(2)(o) 

which may easily be transformed into the more familiar form : 

n(2)(fe2)_n(2)(0) = / dyy{l-y)ln(l + y{l-y)—) (2.16) 

vr Jo ^ rn 

The real beauty of this calculation of this old and familiar result is that the UV divergence 
does not appear until the very last step, the integration over the proper time; and so does 
not disrupt the gauge symmetry of its elements. In contrast, the familiar Feynman graph 
computation in momentum space is so badly divergent that the underlying gauge symmetry 
is lost, and must be reinstated by other means. This was, of course, known to Schwinger, 
who originated proper time calculations in QFT; but it is made clear upon employing the 
elegant and most useful representations of Fradkin. 

3 Radiative Corrections 

In this Section, we display the radiative corrections to the simplest, one closed fermion 
loop, and describe the cancellations which appear even before the corresponding FI is eval- 
uated. These radiative corrections correspond to the action of the linkage operator acting 
upon the omitted A dependence of Section 2, and can be succintly written by inserting the 
terms : 



s 



^ -igo / ds' u Js') A^{x' - u{s')) /go ds' F^„{x' - u{s')) 
e ^ e -'0 tr e .^0 



(3.1) 
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under all of the integrals of (2.6). Eqs (2.7)-(2.9) are still relevant, but before the Jd'^p and 
the Jd[u] are performed, the u dependence of (3.1) must be extracted. For this, (3.1) may 
be rewritten, in the convenient form : 

(ea,-'*!*'"' -') eS. [e^''t.(e'°l''"''^'')Y (3.2) 

where = -^J -D,-- 

The first simplification to be noted was proven in Appendix B of ref.[6], 



n / 90 [ ds'a.F^ 
; ^ tr ( e ■''o 



e'^^trl^e -'o " ) 
and can be trivially generalized to the more relevant statement : 



= tr 1 = 4 



rs rs 

T>. / 90 I ds'a.F. , go / ds' a.F 

e ^ tr e -'o ) = tr ( e -'o 



(e Jo • ) =tr( 



so that the self linkages acting on this OE factor exactly cancel, to all orders in the cou- 
pling. In Feynman graph language, this would correspond to momentum space cancellations 
occurring in every higher order; in the Fradkin representation, one sees them immediately. 



„ -igo f ds'u'.A 
The self linkages e ^ e o 



produce the dependence : 



/ dsi / ds2 u' isi)Dc,^,u{u{si) - u{s2))ul{s2) -igo / ds'u'.A 
e 2 Jo Jo e Jo (3.3) 

leaving the cross linkage operation : 

T.F 



rs rs 

-igo / ds'u'.A / 9o ds' a. 

Jo e tr f e 



(3.4) 



to be evaluated. The A dependence inside the OE can be extracted by writing the latter 
as : 



where x^lu and 5/5x^iu are antisymmetric in ^ and v. In this way, the cross linkages produce 
the term : 

2igl / dsi I ds2u' {si)dxDc,i^u{u{si) - u{.S2)) 



Jo ^ 0X\v{s2) (3.5) 



, i / ds' a^^x^iAs'). 

which is then to act upon the waiting ye Jo j factor. We recall that, in the 

Feynman gauge, Dc,^,,{u{si) - u{s2)) = -^^i^^ [(Au)^ and Au^ = n^(si) -u^{s2). 
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Repeating the argument given following (2.8), one can immediately see that the terms of 
both (3.3) and (3.5) are individually and manifestly gauge invariant. We would like to argue 
that the contributions of the exponential of (3.5) vanishes upon the evaluation of the FI, 
using reasoning very much like that of the Conjecture of ref.^, wherein a rescaling of the 
variable u^{s') — t- \/Qu^{s'), Q real and greater than 1, produces a factor 1/Q multiplying 
the term of that exponential. In [6|, where u{s) + z = 0, a related scaling operation, as 
Q was increased, led to the effective replacement of the term corresponding to (3.5) by a 
gauge dependent exponential factor proportional to In Q appearing outside of its FI. In the 
present gauge invariant case, where z = 0, such a rescaling changes (3.5) to : 



A. 



dsids2 u'^{si) 



Aux 



or to 



— // dsids2U [si] , ■ 12 A- / ^ / i/Q h-X 



(3.6) 



(3.7) 



As Q becomes arbitrarily large, these exponential terms become arbitrarily small, while 
the Q dependence buried in the remainder of the FI is such that the result of the FI, 
with (3.6) or (3.7) replaced by unity, is independent of Q, regardless of the magnitude of 
Q. In other words, such an interchange of limits, taking Q arbitrarily large before the FI 
is evaluated, appears to generate a result independent of (3.5), which would be a most 
convenient simplification. 

We cannot prove this conjecture with any mathematical rigour - and it may be false - 
but for the purposes of this paper it is not necessary. What is immediately clear is that the 
lowest order expansion of (3.6) or (3.7) will contribute nothing to the fourth order estimate 



produces the factor 

2 



of K^ivik)^ because a single bjbx\^ operating upon tr [e. 

order exponential term of (3.3) under the integrals defining KjfJ (k) : 

d'^p 



x^o 

0. K^' is therefore given precisely and exactly by the simple insertion of the Qq 

Z?-{2)/ 



"ds 



-ism 







dsi I ds 







N / d[u] e 2 



d d 



dsa dsb Sgf,{sa) Sgu{sb] 



+ {k^kp S^i^k ) 



■-.ju.U + g) 



u{2hy 

/V I rl\ti\ f> Z, J 

(2^)4 

dsids2 u {si)-u' {s2)Dc{u{si] 



with Dc again in the Feynman gauge. 

All that is now needed is to insert a Fourier representation for Dc{u{si) — u{s2)) '■ 



ffd e^(l-{u{si) - u{s2)) 



le 



dT 



d^q 



-iTq^ + iq.{u{si) - u{s2)) 



(3.9) 



to replace theexp [iq.{u{si)-u{s2))] of (3.9) byexp[z j^l{s')u^{s')] where j^(s') = q^ji[5{s' - 
si)-5{s'-S2)]; to replace the n'(si)-u'(s2) of (3.8) by {d/dsa){d/dsb)[\)^ {5/5j^{sa)){5/5Usb)) 
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and to then perform the gaussian integrals over u and q, before taking the hmits Sa^b — ^ ■si^2- 

rs rs rs rs roo 

The result is a set of integrals over / dsi / ds2 / dsi / ds2 / dr, to be followed by the 

Jo Jo Jo Jo Jo 

f°°ds ■ 2 

final / — e ; and this computation should be vastly simpler than that of the original 

Jo s 

Jost-Luttinger calculation [5j. 

There are two, related situations where linkages operating on OE terms may be expected 
to give a zero contribution to all orders, and these arise in the Dominant Part ( DP ) Model 
of Section 4, and somewhat differently in the Extended DP Model of Section 5. In both 
of these cases we will be interested in a specific limit in which a Au^ variable approaches 
zero, and we will insist that these limits should be taken in a symmetric way, for both 

the odd function dxDJAu) = ott-t — -tk ttt, and the even function dxdn DJAu) = 

i 1 / c- , AuuAua n t 1 • 1- ■ A 

[Ofiu — ^TT-r — 72 r-rj- In such a symmetric limit, as Am — )• 0, we will 



27r2 [(An)2 + ie]2 ^"^-^ [(Au)^ + ie] 
employ : 



dxD^iAu) ^ =0 (3.10) 



and : 



dxd, D.iAu)l^^^ ~ (V - • \ V (A.)^) = (3.11) 

and so replace (3.5) by unity, as well as the cross linkage operation between the OE terms 

4 The DP Model for a single closed fermion loop 

~ (2) 

Return to the KliJ (k) calculated in Section 2, written as : 

Kj^Jik) = ik,K-5,^k')U(^\k') , n(2)(fc2) = 1^ ^^-xsm n(2)(A:2,,) (4.1) 
where : 

n(2) (A:^ s) = ^£dyyil-y)e "^^^'^ " ^) 

We now insert the missing exponential part of this integrand, exp / uDcu] coming 

from all the remaining terms of e'^ — : and retain only those parts which can 

5 AS A A=o 

contribute to subsequent UV divergences; these are the " dominant " parts, which define 
the DP Model, and in the following way. 

As noted in Section 3, the choice of gauge is irrelevant, and we choose the simplest 

Feynman gauge, where the exponential factor — J u'Dcu' becomes : 

-T^ fdsi rds2u{si)-u'{s2)\{u{si)-u{s2)f + iey^ (4.2) 
47r^ Jo Jo L J 
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Since we are concerned in this and the fohowing Section with a particular Au hmit, 

all subsequent OE terms will be discarded. 

It is intuitively clear that the most significant contributions will arise when u{si) is close 
to u{s2), and we therefore expand u{s2) about ,si, writing : 

u^(s2) ^ u,j,{si) - (si - S2)u'i^{si) H 

with u{si) and u'{si) considered as continuous functions. ( The continuity of u{s') is clear 
from its definition, as the integral over the Fradkin 4 velocity v{s'); while the continuity of 
u'{s') = v{s') follows from the physical expectation that the 4 velocity of a particle, real 
or virtual, must be treated as a continuous function of its proper time parameter.) All 
higher derivatives need not be continuous, and there is no obvious way of calculating them 
and their fluctuations; but they should not contribute to the leading divergent structures 
produced by the DP Model. This point is discussed and justified in detail in Appendix B. 

To test the DP Model, in Appendix A we exhibit a completely independent, and simple 
perturbative example, whose log divergence is - to within additive constants - precisely 
the same as that calculated by the DP Model. But this point should be intuitively clear : 
because all u fluctuations are controlled by the gaussian weighting of the Fradkin represen- 
tation, all u fluctuations must satisfy u < y/s; and because all UV divergences arise from 
small s, and therefore from small si — S2, differences which scale as s, we may retain only the 
(si — S2)u' (si) part of the denominator of (4.2), and replace the numerator u'{s2) by n'(si). 
These simple replacements define the DP Model, and effectively permit the extraction of 
the leading divergence structure from under the FI of the u fluctuations. 

The DP Model thus replaces (4.2) by the far simpler quantity : 



47r^ Jo Jo L 



fSl „ .„ .„ —J. 

s' = Si — S2, u''^ = u'^u'^ 



where it is understood that the effective UV cut off e is to be held fixed until the very last 

step of all calculations. If [u'(si)]^ = 0, this integral vanishes; if not, with ao = dol^T^ it 

can be rewritten as : s sx \ 

— — / ds\ \ ds' s''^ + iee{s\) 
TT Jo Jo J 

rsi r- -, _i 

where we denote by e(si) the sign of [u'(s)]^. And since / ds' s' + iee(si) can be 

Jo L -I 

rewritten as / ds' s'+iee(si) = z, the first of the two needed integrations 

Jo L ^ si + iee iee 

can be trivially performed. 

I fs I 1 P _ 

We now insist that the second term immediately above, — / dsi—. — r- = — / dsi e(si), 

le Jo £{s\) le Jo 

must vanish, and for the following several reasons, (a) Intuitively, nothing in the formalism 
distinguishes between positive and negative values of [u'Y-, and therefore the sum over the 
sign of all possible [u']^ fluctuations must vanish, (b) If this integral did not vanish, there 
would appear in all orders of perturbation theory its coefficient which corresponds to a 
quadratic UV divergence; and in all known examples of perturbative, gauge invariant, QED 
calculations, such quadratic UV divergences are absent, (c) One can point to related, if 
indirect, arguments, such as treating the FI of the Fradkin representation as a normalized 
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probability function, and using it to calculate the value of the expected sign of the integral 
of tx'^, which also turns out to be zero. Henceforth, we shall assume this most reasonable 
property, that the integral over the sign of u'"^ vanishes. 

With this understanding, our desired exponential factor is modeled by : 



— / asi s\ + iee{s\) 
Jo 



Since we expect e to fluctuate, as the integral over si proceeds, and since the knowledge of 
£ is only necessary near the lower limit of the s integral, we can rewrite the integral as : 

^r^ = ^ln(^) (4.3) 
which quantity misses a possible — ln(ie(0)) additional term, depending on the sign of 

TT 

e(si — )■ 0). For the moment, we suppress this factor, retaining the obvious In^s/e^ depen- 
dence. 

The exponential of (4.3) is our DP Model of all the radiative corrections ( of this Section ) 
to n(^)(/c^) and yields for the sum of those contributions : 

n(2) (fc^) = 2^ [\ y{l -y)r- i-T'^ e "^^[^o + ^^{1 - y)] (4.4) 
vr Jo Js s Ve/ 

where we have cut off the lower limit of the s integral at e, and will later use the identifica- 
tion £ l/h?. Although (4.4) contains log divergent terms in every order, it is interesting 
to note that these additional radiative corrections conspire to remove the necessity of using 
an £ as the lower limit o the s integral if ao > 0. In any finite order perturbation calcu- 
lation, that lower limit e cut off is absolutely necessary; but since wc have summed over 
an infinite number of graphs, there has occurred a qualitative change in the s integrand, 
such that it was not necessary to introduce this lower limit ( although the need for a cut 
off is simply transferred to large values of s ). There are still divergences, in every order 
( of an ao expansion ), but this effect suggests the possibility that the omitted radiative 
corrections, arising from the linkage of all possible closed fermion loops to those loops under 
consideration above, may tend to suppress the divergences which remain in (4.4). And this 
is indeed the case. 

Restoring the neglected phase factor of (4.3), we calculate Z^^ = 1 + n(0) by rotating 
the contour of the s integration to run along the negative imaginary axis, and so obtain : 



TT/m "0 f°° dr / T \«o/t rao. ... ao . r--rn\\] rAn\ 



and note that the reality of Z3 fixes the choice e(0) = — 1. 
Changing variables, in (4.5), to r = emQX produces : 



n(0) = ^/ dxW'^-ie-^^oa^- ^ / "dxW^ 



or 



3 



n(o)^:, f-o V -lU^ofzrs)^ (4-6) 



1/A2x^? 



3 vmQ^ 
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as long as ao > and (A/mo) 3> 1. Z^^ can now be represented by (4.6); or by an infinite 
sequence of log divergent terms : 



n=l 



(4.7) 



An alternate representation of ^ is obtained by allowing the lower limit of the integral 
of (4.5) to approach 0, which yields : 



n(0) = |2(4)^p. 

Stt V 777,5 / "^0 



' " ' - g-T^aoA-i^ WiL) - r(i + ^) 

3 vrnQ^ TT 



and then : 



7-1 



1 + 



3 Vm^J 



"On 



r(i + ^ 



TT 



(4.8) 



The gauge invariant part of the photon propagator, in this approximation, is given by : 
where : 



-1 



^2 [1 + U{k^)] 



(4.9) 



and n(/c^) is the sum over those proper self energy terms ( which cannot be constructed 
from an iteration over lower order terms ) calculated in DP approximation. The simplest 
renormalization procedure in QED proceeds by adding and substracting 11(0) in the inverse 
of the denominator of (4.9), so that : 

D'^{k) = (A;^)-^ [1 + n(0) + U{k^) - n(o)] ~^ 

and the wave function renormalization constant Z3 identified as the coefficient of the fc^ pole 
of D'^{k) as k"^ — )■ 0, which leads to the familiar identification of Z^ in terms of 1 + 11(0). 
The renormalized propagator is : 



or : 



D.R{k) = Z^'D',{k) 



Dc'Mk) = (k^r' h + Zsiuie) - n(o)) 



(4.10) 



where the combination Z3(n(/c^) — n(0)) defines, in each sequential order, a finite contri- 
bution given in terms of the a constructed from the Z3 contribution of that order, and ckq 
is chosen to have whatever ( large ) value is required so that the renormalized a = 1/137. 
In our DP Model, the calculation of U{P) - n(0) is again immediate, yielding : 



1 



( k^ 
Jo ^ rriQ 



and : 



n(fc2) _ n(o) 



1-6 / dyy{l-y)(l + ^yil 
Jo ^ 



(4.11) 
(4.12) 
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and where all the oq dependence of (4.11) has - in any finite order - been " sequentially 
transformed " into the a dependence of (4.12). 

In standard texts [8J, the finite order passage from (4.11) to (4.12) is emphasized, and 
correctly, as incorrect; but if the conventional, sequential renormalization procedure is to 
be maintained with the inclusion of an infinite number of divergent radiative corrections, 
the situation is quite different. Of course, the conventional procedure of renormalization 
has long been considered somewhat dubious, since it involves an effective perturbative 
expansion in powers of ( the arbitrarily large ) ao- But the procedure does give reasonable, 
renormalized expressions, such as its 27 Mc contribution to the Lamb shift, and the detailed 
calculations [8] showing that the renormalized corrections to this vacuum polarization 
tensor for large (k/m) are proportional to ln{k/m) for both order a and a^, while the 
leading term is proportional to the square of that logarithm. 

But the fact that our DP result for large (k/m) does not agree with the Jost-Luttinger 
perturbative result is irrelevant, for there is no reason to expect the DP model, which deals 
with the extraction of log divergences, to provide the correct limiting values of perturba- 
tive quantities, such as the coefficients of the leading, and finite, a^ln^(A;/m) dependence, 
although such a possibilty may be approximately realized. While our DP Model sums up 
all the possible perturbative divergences to the gauge invariant photon renormalization, it 
misses possibly important additive k/m dependence to those divergent logs, for it is com- 
puting in the region of A; = 0, rather than large k/m. But the DP method does show, as 
in Section 5 for the A: = computation of Z^^, that it is sensible to consider this gauge 
invariant sector of QED as a finite QFT; and this, conceptually, is a new and most satisfying 
result. 

Finally, to understand what it is that our DP Model achieves, we ask and answer 
the following question. How can one understand the connection between the vanishing 
of a propagator's denominator in this configuration space, functional formulation and the 
appearance of log divergences in momentum space ? Simply by taking the Fourier transform 
of Dc{u{si) — u{s2)) and asking why does the corresponding 



converge when combined with another momentum-space propagator? It is not the insertion 
of a denominator factor proportional to [{q — fc)^ + M^]~^ coming from another propa- 
gator, bosonic or fermionic, for that alone generates a log divergence. Rather, it is the 
exp[ik.{u{si) — u{s2))] term - which of course is eventually evaluated in and by the Jd[u] 
FI - which provides enough oscillations and cancellations to yield a finite result. But when 
u{si) — u{s2) vanishes, as we have suggested above, that log divergence will appear and will 
enter all relevant parts of the computations of that order. We find it far simpler to remain 
in functional configuration space, where the cancellations of the OE terms are easily visi- 
ble, then to convert to the conventional Feynman graph analysis in momentum space. The 
divergent logs will reappear, but the cancellations obvious in functional configuration space 
will there require much tedious calculation to obtain. It is far more efficient to identify and 
extract in their functional source those log divergences, than to perform the FIs, convert 
to momentum space, distribute those divergences in the conventional parts of a Feynman 
graph, and then attempt, order by order, to understand and remove them. 
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5 The extended DP Model, and the finiteness of Z-x 



We return to (2.3) and the exact K^^y, from which follows the complete ^ . The first 
task to perform is to argue that the : 



8L 5L 
5A^{x) 5A^{y) 



/<S> 



A=0 



terms cannot contribute to Z^^ , and for this we return to the Fradkin representation for 
L[A], here simplified by the neglect of all OE terms, as discussed in the previous Section. 
Together with the Proof of Appendix B, this means that all OE dependence of every L[A] 
does not contribute to the radiative corrections calculated in the DP Models. 

SL SL 

Consider first the product , , , . , which contains under its separate integrals 

oAi^[x) dAj^[y) 

the terms relevant to this discussion: 

dsi Un{si)S{x' — u{si) — x) jd!^x" 5(u{s)) J dsi Ui,{si)d{x" — u{si) — y) 



X exp 



rs rs 

-igo ds u'^{s') Aa{x' - u{s)) - igo ds Upis) Ap{x" - u{s)) 



Neglecting for the moment cross linkages to the exp(L[^]) term, we then have : 



5A^,{x) 



6A,{y) 



A=0 



/rs rs 
d'^x" 6{u{s))j dsi J dsiu'ij^{si)ul{si) 6{x' — u{si) — x) S{x" — u{si) — y) 



X e 



J Jo Q I J Jo 

igo / ^^^1 / ds2 u'{si).u'{s2) Dc{x' — x" + u{s2) — u{si)) 
X e .^0 Jo (5.1) 

With the aid of the delta functions of (5.1), its last line may be rewritten as : 

rs rs 

igo / dsi ds2u{si).u{s2)Dc{x-y + A{u,u)) 

(5.2) 



JO 



and the Jd^x'Jd'^x" performed so that (5.1) reduces to : 



d{u{s))S{u{s))J^dsiu'^{si) j^dsi <(si) 5(An)5(Au) 



X e 



igl / dsij ds2u'{si).u'{s2) Dc{x - y + u{s2) - u{si)) 
Jo Jo 



(5.3) 



where A{u,u)) = u{si) — u{s2) — u{si) + u{s2), and S{Au) and S{Au) refer to the self 
linkage exponentials of (5.1). 
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The renormalized charge in QED is conventionahy defined by evaluating Kfj_y{k) at 
= 0, corresponding to the definition of physicaUy measured charge at large distances, 
specifically at distances large compared to the Compton wavelength m~^; the radiative 
corrections occur at distances less than m~^, and in the context of these calculations, this 
corresponds to evaluating Kfj_u{x — y) at separations x — y » m^^. But all of the x — y 
dependence of (5.3) lies in the corresponding to its cross-linked exponential factor, while 
the u and u quantities, from their definitions scale as -y/s ~ \/l ~ . It is then clear 
that in the limit of x — y >> m~^, the An and An dependence of Dc{x — y + A{u,u)) is 
effectively suppressed, and in this limit the factors Jq dsiu' (si) and ds2v! {S2) = 0, so 
that this exponential factor completely disappears. In a similar way, after the x' and x" 
integrations have been performed, the entire contribution of (5.3) is itself proportional to 
similar factors, dsiu'^{si) dsiu'y{si) which also vanish. 

Insertion of the cross linkages between (5.1) and the A dependnce of exp(L[A]) does not 
change this situation, for the limit of large x — y'm. the cross linked terms can be understood 
as the separate limits of x — )• 00 and y — )• —00, so that these cross linkages also vanish. The 
computation then reduces to that of (5.3), as the self-linkages of exp(L[74]) are cancelled by 
the definition of < 5 >. The result is that the entire quantity : 



6A^{x) 5A^{y) 



A=0 



does not contribute to Z3. 

The first term of (2.3) is the relevant quantity, and we first consider the structure of 
the factor c^C^L It will be convenient to express the latter in terms of the functional 
cluster expansion [9j : 

00 -| 

e^Aei[^]=exp[^lQ„[^]' 
L-' — ' n J 



n=l 



where Qn[A] = e^A(L[A]Y , that is : Qi[A] = e'^^L[A] = L[A], Q2[A] = 

V / connected 

(e-^L)(e'^A — l){e'^L), etc. The reason for choosing this expansion is that the Qn, 
n > 1, can be estimated to yield smaller values than does Qi (as well as being far more 
difficult to calculate); their divergence structures are similar to that of Qi, but, as discussed 
in Appendix C, they play a smaller role in the overall calculation. Qi alone is sufficient to 
remove all the perturbative log divergences of . 

Rather than repeat all the details of every equation in the next few paragraphs, we 
shall simply present the added features that arise from the cross linkages between the 

e — , , ^ , — —r of the previous Section and the e-^e'^~e^of the present discussion. 
6Af,{x)SA^{y) 

Note that all the self linkages of e ^ e 

L[A] 

/< S* > simply disappear from the final result 
by virtue of the definition of < S" >. We are therefore interested in : 



e^A il ]^Va,L[A] 

6Af,{x)5A^{y) 



(5.4) 



which, continuing to use the variables u„(s') for the FI of e-^ -— — — -— corresponds 

dAf,{x)dA,y{y) 
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to the insertion under the latter's FI the quantity : 



exp 



-50 /W / d^w Dc{x' - u{s') -w)—^ 

Jo J OAa 



A-^O 



(5.5) 



where we again hold to the Feynman gauge. 

But the operation of (5.5) is just a translation operator, and regardless of what it acts 
upon, has the effect of shifting the A dependence of that function - in this case L[A] - by 
the quantity : 



exp 



iglTds' <(s') I'dt' v'^{t') D,{x' - x" - u{s') + v{t')) 
Jo Jo 



(5.6) 



appearing under the FI of L[A], with x" and Va{t') the variables of that functional. In 
addition to the factor of (5.6), there appear under the L FI the self linkages of amount : 



exp 



JJl^dtidt2 DMti)-At2))v'{h).v'{t2) 



which are independent of x' and x". In effect, what the cross linkages have achieved is to 
insert s dependence under and mixed with the t integrals of L; and subsequent integration 
of that t dependence generates an exponential s dependence which will have a damping 
effect on all the divergent expressions of Section 4. This simple observation is at the heart 
of the mechanism for obtaining a finite charge renormalization. 

The self linkages of L may be read off from those calculated in Section 4, eq.(4.3) and 
its subsequent discussion : 



exp 



Jjytidt2 D^{v{ti] 



V{t2))v\ti).v'{t2 



The cross linkage term of (5.6) can be evaluated in a similar manner, with the realization 
that only the continuous parts of the functions u{s') and v{t') are relevant - as noted 
in Appendix B - and in essence they are very similar functions, differing mainly in their 
physical place in the calculation. They represent the continuous parts of the fluctuations 
defined by the same functionals; and it is reasonable to ask when they can interact and 
combine directly with each other, for when u{s') ~ v{t') one has the beginning of an 
incipient divergence. A modification of the previous DP model is now defined by expanding 
v{t') ~ v{s') + {t' — s')v'{s') + • • •, for the case when these functions are essentially the 
same: when u{s') c± v{t') and u'{s') c± v'if), so that the denominator of the cross linkages 
becomes : 



{{x' -x") + {t' - s')v'{s')f + ie 



(5.7) 



or the same form with v replaced by u. Since the functional integrations sum over all 
possible (continuous) forms, which fluctuations are defined in exactly the same way, there 
should be a strong possibility of such overlaps. The essence of this Extended DP Model is 
that u{s') and v{t') are treated in a completely symmetric manner; and this requirement of 
symmetry turns out to be a guarantee of simplicity of the forms that follow. This intuitive 
assumption defines the Extended DP Model, in which the difference of two, equivalent. 
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continuous fluctuations of identical functionals has the same possibility of overlap as in the 
DP Model. We emphasize that we cannot prove the validity of this assumption; but it is 
most certainly intuitive; and it forms the basis of the cancellations of divergent logarithms 
that are about to occur. 

However, even the EDP Model cannot guarantee the vanishing of (5.7), for another 
difference, x' — x" , need not be small, and in most cases is not. What is the consequence when 
x! —x" is large ? Quite independently of the EDP Model, when x' —x" » m^^, the Conipton 
wavelength of the charged fermion traveling about the loop, that difference completely 
dominates u{s') — v{t'), because each of the latter quantities scale as s/s ^ \/t ^ m~^, and 
hence the u — v difference in that case is irrelevant. But then, as repeatedly emphasized, 

the s' and tf integrals of (5.6) vanish, / ds u^{s) / dt' v'^{t') = and remove the entire 

JO Jo 

cross-linkage term from consideration. When does this not happen ? Only when x' — x" is 
restricted to values on the order of, or less than u{s') — v{t'). But we are interested in small 
differences, where s' and t' tend to the order of e, and where subsequently e — )• 0. How can 
this be arranged ? 

We shall here assume that the only contribution to the cross linkage integral comes 
when \x' — x"\ ~ where ^ ~ 0{l). Since this is an idealization, one must expect 

fluctuations about this condition, such that ^ will turn out to be somewhat less than unity. 
For conceptual simplicity, choose the point x' as the origin of the x" coordinates, and 
consider a (Euclidean) 4-sphere of radius ^-^/e For any point within this sphere, \x' — x''\ 
effectively disappear from (5.7), and we can apply the EDP Model; this means that the 
only non-zero values of the x" integrals is given by : 

Can (5.8) produce a non zero result ? Yes, because L is itself proportional to the fac- 
tors -^tr[l] jd'^x" j^je-^^"^^ Njd[v] j v{2h)-^v]6^^'^ (v{t)) and, just as for 

the DP Model, the evaluation oiNjd[v\ yields : N j d[v\ ex.p[^ j v{2h)~'^v]6'^'^^ (v{t)) = 
7r2 1 

—i^-, — ttj so that the entire set of exponential integrals reduces to : 
i2 (27r]4' ^ ^ 

^(?) =4)'"' f^e-iirnl^3inp/2 (J)^(i)^^ (5.9) 

Jttp/2 fi^P 



In (5.9), one factor of e ^' [^-j arises from the self linkages of L, while the cross linkages 
generate the remaining factors. We have maintained strict s, t symmetry by writing : 



t ( rs rt r ft(e'-t') 9{t' - s' 



ds' / dt' , , ^ , ,,\ ; . = ds' dt' 
{u{s') -v{t'))^ + i£ 







/•s r 1 11/"* 
= / ds' ——- - + / dt' 

Jo ll££{s') S' + I££{s')j Jo 



t' + ieeis'))^ (f - s' + iee{t'))^ 
1 1 



1 1 

4- / H+' 

i££{t') t' + i££{t') 

Again, the integrals over the £ factors vanish, and the result is simply : 



\i££{0)J Vi££(0)/ 
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which is properly symmetric in s and t. Were that symmetry not preserved, the results 
would lead to far more complicated forms, requiring detailed numerical integrations in 
order to verify the expectations of a finite and an a close to 1/137. In contrast, the 
symmetric Extended DP Model adopted here leads to results obtainable in closed form, and 
to the immediate verification of our expectations. 

With the inclusion of these cross linkages, the DP integral representation for Z^^ is 
changed to : 

Note that we are keeping to the conventional perturbative form (although in configuration, 
rather than momentum space) of cutting off all proper time integrals with a lower limit of 
e, which will shortly be set equal to zero. Without the cross linkage factors that produce 
r, (5.10) has the same divergences as does (4.5); but an entirely new situation now arises 
with the insertion of (5.9) into (5.10). Proper time contours need not be rotated; all that is 
needed is the simple change of variables: s/e = x, t/e = y, so that (5.10) may be rewritten 
as : 

^3-1 = 1 + (|) e'^'P/'^ j^dxxP-^e-^^'^l^e^^^) (5.11) 

where : 

r(a;) = i(i)'e3z7rp/2^p j^<^ ^-ieyml y2p 

and one notes that, for convergence as y ^ 00, one must have p <1. For p > 0, T{x) can 
act as a damping or oscillating factor that provides convergence for the x integral as — t- 00; 
and for both the x and y integrals, we may now safely let £ ^ 0, and the divergences have 
disappeared. 

As this program is carried out, one notes the independence of Z^ on mo, or on m. 
This property is not at all clear from perturbation theory, where sequential renormalization 
of mass and charge, along with simultaneous changes in Z\^2 must appear. In fact, Z3, 
and therefore a, are independent of the charged particle's mass, in agreement with the 
experimental fact that all charged fermions obeying QED (but not simultaneously QCD) 
have the same electric charge. 

What remains is to insure that Z^^ is real, and hence the condition Im Z^^ = specifies 
a relation that p = uq/tt must satisfy. In principle, this is true; and if that condition leads 
to a single allowed value of ao, and hence of a, one will have solved an old and deep question 
in Physics. In our calculation, however, there appears the parameter ^, as a measure of 
the difficulty and uncertainty of extracting the divergent character of (5.7). And with the 
neglect of higher terms of the cluster expansion, there is no guarantee that a single value 
of ao will emerge. Rather, with the cluster approximation already made, and with those 
approximations we are about to make, we are gratified to find a range of p values, Q <p < 1, 
within which we can choose ao such that a ~ 1/137. 

For p < 1, the integral defining T{x) is readily obtained, and the expression for Z^^ 
now reads : 

= 1 + (|) e '^'^pI'^ j^dx xP-^ e -^^^ (5.12) 
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where the q parameter is given by : 

.^s 4 e3«7rp/2 

Note that convergence of the x integral is obtained for Re(g) > 0, which corresponds to 
p < 2/3, and that a necessary but not sufficient condition for the removal of ImZ^ is that 
p > 1/3. 

The X integral is completely trivial because its integrand is a perfect differential, and 
one finally obtains : 

Z^' = l + l{l-p){^yUe-'-Pe-^ (5.13) 

Since we expect ^ to be somewhat less than 1, the ration (^/2)^ should be very small, and 
a good approximation to (5.13) is then : 

= l + l{l-p) (i)"' (sin(7rp) +icos(7rp)) (5.14) 

and it is clear that p = 1/2 insures that the Z^'^ of (5.14) is real, while the choice ^ ~ 0.397 
leads to an a = npZ^ of approximately 1/137. Prom this solution, one calculates (C/2)'^ ~ 
0.00155, so that any correction to these parameters obtained by the use of (5.13) rather 
than (5.14) cannot differ from the above values of p and by more than a few parts per 
thousand. 

6 Summary 

The thrust of this paper has been to argue, by summing the " naturally divergent " 
terms of all relevant radiative corrections, that charge renormalization in QED is finite. We 
do not claim to have given a mathematically rigorous proof of that statement, but rather 
an intuitive statement, based on the functional structure of QED. We have argued that 
our extraction of logarithmically divergent terms corresponds to those found in lower order 
radiative corrections using Feynman graph techniques; and we believe that in momentum 
space, graphical techniques become impossible, and therefore irrelevant, in any attempt to 
include all, or almost all, radiative corrections of arbitrarily high order. 

The functional techniques we use are based upon a convenient rearrangement of the 
Schwinger/Symanzik functional solution for the generating functional of QED, together 
with a slight rearrangement of Fradkin's most useful functional representation for Green's 
functions Gc[A] and, in particular, for L[A], the log of the so called " fermion determinant ". 
L[A] contains the basic, gauge invariance of the photon propagator, and that structure is 
here realized by means of most convenient linkage operations. 

If a perturbative expansion is desired, this functional approach will exactly reproduce the 
conventional Feynman graphs, but it has the great advantage of working in configuration 
rather than momentum space, and one can take advantage of cancellations which occur 
there before any computation is required, but which are achieved in momentum space only 
after painful and tedious manipulations [Tj. And, it should be noted, that frequently, as is 
the case for the lowest order radiative corrections to the photon propagator, singularities 
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of the Feynman integrals can mask the symmetry structure of the theory, an unpleasant 
attribute of that method of calculation which is quite absent from functional methods built 
around proper time representation^. 

Within this formalism, and the intuition we have used to obtain the results described 
above, we have derived a pair of integrals in (5.11), with e set equal to zero, which should 
express the finite character of charge renormalization. And we have illustrated a possible 
solution of those equations with a simplified model which generates a specific value of uq 
that leads to an a ~ 1/137. What this model solution does not do is to obtain a single, 
precise and necessary value of uq, and so determine a. As noted in the text, because of 
the imprecision in extracting the divergences from the cross linked, closed fermion loop 
functionals, as expressed by the parameter ^, this goal may be elusive; and we will have to 
be content with choosing an uq which does reproduce the experimental a. 

Previous attempts have been made some decades ago [1] to simultaneously display a 
cancellation of divergences, which might lead to a value of a close to its renormalized value. 
These were noble efforts, especially in the context of Feynman graphs; and one can now 
see why they were unsuccessful, for the crucial aspect of including an infinite number of 
closed fermion loops, each containing all possible photonic "dressing" in a manifestly gauge 
invariant way, could not be done. It is gauge invariance, built into the Fradkin representation 
for L\A\, along with the use of proper time techniques which preserve that invariance, and 
allows one to identify and extract divergences. 

It should be noted also that the assumptions made in these papers, which center about 
a single log divergence for Z^^ persisting in higher perturbative orders, together with the 
special choice of a zero bare fermion mass, are quite different from what we observe. Our 
Z^^ is finite and independent of mass. The use of a finite number of terms in a Bethe- 
Salpeter kernel (cf ref[10]. Section 1) misses an infinite class of Feynman graphs, for a 
new class of graphs appears in each higher order; and in each missing class there are an 
infinite number of Feynman graphs for higher and higher orders. In contrast, we make no 
assumptions about the nature of perturbative divergences in Z^"^; we include them as we 
see them, and we are able to see and include them because we use a formulation where their 
appearance is obvious, and where strict gauge invariance helps us to extract them, and sum 
them to all orders, literally. 

There are other, relevant questions which immediately come to mind, such as why one 
should expect the approximate value of a to be close to 1/137, when no mention has been 
made of electroweak symmetry, in which radiative corrections involving the weak interac- 
tions should play a role. Were we calculating a process at significant momentum transfers, 
or at energies on the order of the W mass, for example, then our simple restriction to 
electron and photon QED would be insufficient. But the Z^ computation involves only 
radiative corrections evaluated at A; = 0, a most non-relativistic limit; and even though 
UV divergences may appear in every perturbative term, when summed, the result is inde- 
pendent of mass, as it becomes finite. Because the weak interactions are far weaker than 
electromagnetic interactions, we do expect corrections to our approximate calculations of a 
to be relatively small. 

^Functional representations provide a simple realization of the " Last Rule ", often stated and rarely 
understood, for writing the total of all Feynman graphs of a given order: " Sum over all topologically 
distinct graphs ". 
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It might also be noted that, from general principles enunciated long ago by Schwinger, 
one expects the renormalized a to be smaller than the unrenormalized aoj and since our 
result for the latter is appproximately 7r/2, it is a happy circumstance that the choice of 
^ ~ .4 - which, as expected, is on the order of but somewhat less than 1 - does produce an 
a close to its measured value. 

Another immediate question concerns the possibility of employing the present methods 
to evaluate Z2; can the analysis of the previous sections suggest, for any choice of relativistic 
gauge, that is finite ? Unfortunately, this does not seem to be possible, and the 
reason is that the conditions stated in section 5 for the validity of the EDP Model are 
no longer satisfied. Instead of the properties u(s) = n(0) = f (t) = u(0) of the photon 
calculation, which allowed one to expect the similarity of the functions ti(s') and w(t'), for 
the electron propagator one has n(0) = v{i) = v{0) = 0, but u{s) = —z = —{x — y), 
the variable conjugate to the electron's momentum. And as one goes to the mass-shell in 
p^, one expects z'^ to become large, so that there will be a considerable difference in the 
(continuous) variations of u{s') and v{t'). 

For this reason, the intuition of the EDP Model disappears, and there remains no obvious 
mechanism to cancel the divergencies, as was done by the function T(x) for Z^^. On the 
basis of this argument, we see no alternative to the conclusion expressed long ago by Kallen 
|10j . who conjectured that at least one of the renormalization constants of QED is divergent. 
From our perspective, it is the gauge dependent, unmeasurable which diverges, while 
the gauge independent Z^^ is finite. 

Perhaps the best way to end this paper is by emphasizing that while the functional tools 
used are surely powerful, and appropriate for this problem, the identification and extraction 
of divergences we have used has been intuitive. Nevertheless, the obvious advantages of this 
functional approach to the calculation of radiative corrections seems clear, and, we believe, 
deserves strong emphasis. 
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Appendix A. Equivalence Example of the DP Model 

To demonstrate the accuracy of the DP Model in a quite different setting, consider the 
calculation of the simplest fermion self energy graph, where - to make connection with the 
boson nature of the L[A] calculations - we simplify by passing to a bosonized version of 
relevant fermion propagator at the same stage of each calculation. This example shows that 
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the log divergences of both calculations are exactly the same. 

The order contribution to the fermion propagator is given, in functional notation, 

by: 



^2 '" 



= {-i){igY / D,{u - w) u I A)7^ G,{u, w \A)j^ G^iw, y \A) ^ (A.2) 

where we again adopt the simplest Feynman gauge Dc^^u — ^fiu-^c 

For the Feynman graph calculation, we replace both of the external Gc by their free 
field limit Sc, and, to have a clear correspondence with the DP calculation, imagine that 
the central Gc is replaced by a bosonized version, < u [[m^ + (9 — i5f^)2]~^| w >, whose 
free particle limit is simply <u \[m? + S^]"^! w>. Taking Fourier transforms, one finds the 
corresponding contribution to (A.2) : 

^9^ /■j4,. 1 1 



Sc{p) If^ d k 2 '^c(p) 



(27r)4 ^'^' '^7 k"^ -ie {p-kY + m? 



or, with Sc{p) = — ^—j , V 7m7m = 4 : 
ni + iff ^ 



d'k ,^ (A.3) 



47r^ Vp^ ) J k'^ — ie {p — ky — ie 

where, for simple comparison with the DP Model, we choose m = 0. 
Using the simplest regularization possible : 

1 1 1 f^\, 21 



dl 7775 ^7777 , A +00 



fc2 k^+ A2 7o {k^ + 

1 /■! _„ 

and the Feynman combinatoric —77 = 2 / dx{\ — x)[a + x{h — a)] , the integral of (A.3) 

a^o Jo 

rl ^2 

is easily shown to be : ivr^ / dx\i\.\ — ^ ). 

Jo ^xp^J 

The result of this computation is then : 

ff2 1 fl .A2 



/ dxhii (A 4) 

47r2 p2 _/q V xp^ > 

We now calculate the equivalent quantity using the DP Model of the fermion propagator. 
We recall the expression for the electron propagator [3j : 

/■oo -2 [■ — u(2h)~^u / 
Gc{x,y\A) = iJ dse-^^"^ Njd[u]e'^J i^m-j- 

X 5^"') {u{s) + x — y^ e 



6u'{s) 

ig I ds' u^{s') A^{y - u{s')) f g j ds a^^ Ff,y{y - u{s)y 
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We suppress all spinorial terms, set m = and z = x — y . The order contribution to 
the fermion propagator is then : 



i j dsN d[u\ e 



2 47r2 



(5^^) (u{s) + z) 
/ dsi / ds2u'{si)-u'{s2) (■u(si) - ■u(s2))^ -f- 

Jo Jo L 



(A.5) 



Using the notation and operations as defined in the text, the DP Model replaces the second 
line of (A.5) by : 



si + iee iee. 



We again argue that / dsi £ = 0, and replace the integral of (A. 6) by : 
Jo 

Jiee{o)Si \iee{0)J 



(A.6) 



(A.7) 



The first line of (A.5) is then immediately given by : i ds 

Jo 

combination with (A.7), its Fourier transform then yields : 



— i I dse~^^P In , 

47r2 Vn \iee{0) 



d^n ■ ■ 2 



(2vr) 



(A.8) 



Rotating the s integration contour so that s — )• — ir, and with £ = A ^, e(0) = —1, x = Tp^, 
this becomes : 



and one sees that the log divergent terms of (A. 9) and (A. 4) are exactly the same. 



(A.9) 



Appendix B 

We can attempt a heuristic justification of the DP Model in the following way. Because 
the Pradkin v{s'), < s' < s, refers to the four velocity of a virtual particle ( e.g., the 
fermion moving in a closed loop of the radiative correction corresponding to the simplest 
vacuum bubble ), wc demand that v{s'), and therefore u'(,s'), be a continuous function of 
its proper time, s'. This is a physical restriction on the class of functions u{s') allowed; 
whether the particle is real or virtual, its four velocity should be and will be assumed to be 
coiitimious. 

From its definition, u{s') = ds"v(s"), u{s') itself must be a continuous function of 
s'; and hence we have restricted consideration to the class of functions u{s') which are 
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continuous and have a continuous first derivative. But no statement can be made about 
higher derivatives, which must be expected to be discontinuous. One way to describe such 
functions is to imagine that they begin hfe as continuous, with continuous second and higher 
derivatives proportional to finite constant Pn- Then, as the functional integration proceeds, 
imagine that certain of these parameters p„ are changed, in a random way, to have values 
which approach ±00; and this corresponds to the introduction of discontinuous higher order 
derivatives. What will then be the effect of such fluctuations on the relevant functions : 



rS /•» p 

/ dsi / ds2u'{si)-u'{s2) {u{si) - u{s2))'^ + ie 



i 

2(27r)4 Jo Jo 

of our functional integrand ? 

Consider first the denominator {u{si) — u{s2)) + is . Before the fluctuations are 
allowed to destroy the continuity of the second and higher derivatives, the main contribution 
of this denominator will, clearly, come from the expansion we have used, replacing {u{si) — 
u{s2)) by u'{si){si — S2), and neglecting terms such as (si — S2)^u"{si)/2. If the pn are 
then allowed to fluctuate such that any higher derivatives approach ±00, the denominator 
will then become infinite, and so conveniently removes itself from the FI. 

Now consider the u'{si) ■ u'{s2) numerator term, which we have replaced by ti'^(si). 
Again imagine that all derivatives begin life as continuous, and consider the first correction 
to this approximation, (s2 — si)u'{si) ■ u"{si). If u"{si) is finite, in the neighborhood 
S2 — si, where the denominator is about to vanish, then this term generates a relatively 
small correction to our approximation; but if and when u"(si) becomes infinite, there will 
be no contribution to the FI because the denominator fuctuations will involve u"{si)^, 
which is more divergent than the numerator it'(si) • u"{si), and hence, both numerator and 
denominator contributions from such a discontinuous second derivative are removed from 
the FI. 

This argument can be repeated for every higher order derivative, and for the sum of 
all such higher order derivatives; and the result is the DP Model defined in the text. This 
argument is heuristic, a physicist's argument; it is intuitive, rather than mathematically 
rigorous. To justify this intuition, we point to the example of Appendix A, wherein the 
log divergence calculated from the DP Model yields exactly the same result as that ob- 
tained from the corresponding Feynman graph integral. In this very real sense, a divergent 
Feynman graph in momentum space, calculated from continuous, if overlapping, momen- 
tum space integrands, may be thought of as equivalent to the " continuous " elements of a 
functional integrand in the DP Model. 



Appendix C 

There are two distinct reasons why the connected Qn, n > 1, may be expected to 
generate a contribution to T{x) which is considerably less than that of Qi. The first is 
because of the reduced probability of finding overlaps of three or more coordinate systems 
whose origins must lie within a distance ( in Euclidean 4 dimensions ) of ^^/£ from each 
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other, one of those origins being that of the x' of 



6'1 



. The probabihty of such 



overlaps can be thought of as proportional to the overlapping volume divided by the total 
volumes of all three or more spheres each of radius ^v^, and that ratio is a small number, 
becoming smaller as n increases. In mathematical terms, a non zero intersection of the 
supports of three or more independent distributions such as DcDc ■ ■ ■ Dc is much less than 
the corresponding quantity for the case of two such independent distributions. 

The second reason is more closely tied to the computations of Section 5, and has as 
its origin the nature of the " connectedness " requirement, which can be stated in the 
following way. Connected linkages require that there shall be at least one linkage between 
the connected parties, e.g., for Q2 : 



Lie 



-Da 



1]L 



(C.l) 



with the result that the linkages between the factors of exp 



-igo fdt'v'{t')-A{x"-v{t')) 

Jo 



from one L, and the factor exp 
appear in the form : 

rt 



-igo / dr'w'{r')-A{x"'-w{r')) 
Jo 



from the other L, will 



exp 



iQo [ dt' fdr' w'{r')-v'{t') D^{x" - x'" + w{r') - v{t')) 
Jo Jo 



which can be rewritten as : 



J^dX-^exp i\gl J J w'-v'Dc 



or as 



Wo 



I I r-i 



dX exp 



w 



■V Dr 



(C.2) 



Assume that the overlaps have taken place and consider the multiplicative term of (C.2). 
It does not invovlve yi^2 dependence raised to a power, but rather the logarithm of that 
dependence, which under the same variable changes as used in Section 5, will convert to 
X and ui^2 inside logs. From our over simplified model of Section 5, one sees that the 
significant x value of the final Z^^ integral is just barely larger than 1, and hence the ln(x) 
terms of this log dependence will not contribute significantly. The remaining ui^2 factors 
must be evaluated within the / dui^2 of each L. For 1 < ui^2 < x, we get another ln(x), but 
for < 2 < 1, there will be non zero contributions, integrable quantities of 0(1), relative 
to the results of the / dui^2 integrals without such terms. 

However, these log terms multiply an integral which is essentially an average over values 
of < A < 1, over the forms {q)'^^, where q takes on the different values which can be read 



off from the text, e.g. (yi)^f ) ' . Since p has been here replaced by Xp, for 

^yl-y2-^ ^ 

small A the contributions essentially disappear; and only when A ~ 1 is there a significant 
value to the integral. This represents an effective decrease of the effectiveness of the coupling 

between the two Ls; and when multiplied by the small log terms discussed above, in addition 
to the small overlap factors, it seems clear that the connected terms cannot significantly 
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add to T{x), and therefore cannot significantly change the result of the Qi computation. 
Without a detailed and rigourous analysis, one cannot be absolutely sure; but this is our 
(intuitive) belief. 
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